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Abstract. We describe aspects of the role that diquark correlations play in under- 
standing baryon structure and interactions. The significance of diquarks in that appli- 
cation motivates a study of the possibility that dense hadronic matter may exhibit di- 
quark condensation; i.e., quark-quark pairing promoted by a quark chemical potential. 
A Gorkov-Nambu-like gap equation is introduced for QCD and analysed for 2-colour 
QCD (QC2D) and, in two qualitatively different truncations, for QCD itself. Among 
other interesting features, we illustrate that QC2D with massive fermions undergoes 
a second-order transition to a superfluid phase when the chemical potential exceeds 
ra^/2. In the QCD application we illustrate that the a := — (gq) 1//3 7^ phase, which 
determines the properties of the mass spectrum at zero temperature and chemical po- 
tential, is unstable with respect to the superfluid phase when the chemical potential 
exceeds m 2 a, and that at this point the diquark gap is large, ~ a/2. The superfluid 
phase survives to temperatures greater than that expected in the core of compact stars. 
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1 Diquarks 

A diquark is a bosonic quark-quark correlation, which is necessarily coloured in 
all but 2-colour QCD (QC2D). Therefore, in the presence of colour-confinement, 
diquarks cannot be directly observed in a N c > 3 colour gauge theory's spec- 
trum. Nevertheless evidence is accumulating that suggests confined diquark cor- 
relations play an important role in hadronic spectroscopy and interactions. 

The first discussion of diquark correlations in literature addressing the strong 
interaction is almost coincident with that of quarks themselves It was 

quickly realised that both Lorentz scalar and vector diquarks, at least, are im- 
portant for baryon spectroscopy ]3| and, from a consideration of baryon magnetic 
moments [|J, that the diquark correlations are not pointlike. This latter point 
is still often overlooked, although with decreasing frequency and now certainly 
without the imputation that it is a realistic simplification. 

The motivation for considering diquarks in the constituent-quark model is 
that treating baryons directly as a three-body problem poses significant chal- 
lenges in anything other than a mean-field approach. The task is much simplified 
if two of the constituents can be replaced by a single degree of freedom. However, 
an obvious question is whether there is any sense in which that replacement is 
more than just an expedient; i.e., a sense in which it captures some important 
aspect of QCD's dynamics? The answer is "yes" and we now turn to explaining 
that. 

A significant step toward a description of baryons in quantum field theory can 
be identified in the realisation M that a large class of field theoretical models of 
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the strong interaction admit the construction of a meson-diquark auxiliary-field 
effective action and thereby a description of baryons as loosely-bound quark- 
diquark composites. This is the class of theories with a chiral symmetry pre- 
serving four-fermion interaction, which includes, e.g., the Nambu-Jona-Lasinio 
model H and the Global Colour Model |t]], that have been widely used in 
analysing low energy strong interaction phenomena. 

The picture of a baryon as loosely bound quark-diquark composite can also 
be reached via a direct analysis of the bound state contributions to the three 
quark scattering matrix. The associated Schwinger function (Euclidean Green 
function) is just that quantity whose large Euclidean-time behaviour yields a 
baryon's mass in numerical simulations of lattice-QCD. Considering the colour 
structure of this Schwinger function, we focus on the Clebsch-Gordon series for 
quarks in the fundamental representation of SU C (3): 

3 C <g> 3 C (g> 3 C = (3 C © 6 C ) <g> 3 C = l c © 8^, © 8 C © 10 c , (1) 

from which it is clear that a colour singlet 3-quark contribution is only possible 
when two of the quarks are combined to transform according to the antitriplet, 
3 C , representation. This is the representation under which antiquarks transform. 

Single gluon exchange is repulsive in the 6 C channel but attractive in the 
3 C channel. It is this feature that underpins the existence of the meson-diquark 
bosonisation referred to above. One way to see that is to realise that the auxiliary 
field effective action obtained for any element of the class of four-fermion interac- 
tion models provides a Lagrangian realisation of the rainbow-ladder truncation 
of the Dyson-Schwinger equations (DSEs) The rainbow-ladder truncation 
has been widely and successfully employed in the study of meson spectroscopy 
and interactions, see, e.g., Refs. P, ^0| , pl|^2| , and nonpointlike colour-antitriplet 
diquark bound states exist in this truncation of the quark-quark Bethe-Salpeter 
equation (BSE) Hence they provide a real degree of freedom to be used in 
the bosonisation. 

At first sight the existence of colour-antitriplet diquark bound states in these 
models, and in the rainbow-ladder truncation, appears to be a problem because 
such states are not observed in the QCD spectrum. However, as demonstrated 
in Refs. [|l4j[l5|], this apparent lack of confinement is primarily an artefact of the 
rainbow-ladder truncation. Higher order terms in the quark-quark scattering 
kernel, the crossed-box and vertex corrections, whose analogue in the quark- 
antiquark channel do not much affect many of the colour singlet meson chan- 
nels, act to ensure that the quark-quark scattering matrix does not exhibit the 
singularities that correspond to asymptotic (unconfined) diquark bound states. 

Nevertheless, such studies with improved kernels, which do not produce 
diquark bound states, do support a physical interpretation of the "spurious" 
rainbow-ladder diquark masses. Denoting the mass in a given diquark channel 
(scalar, pseudovector, etc.) by m qq , then £ qq := l/m qq represents the range over 
which a true diquark correlation in this channel can persist inside a baryon. In 
this sense they are "pseudo-particle" masses that can be used to estimate which 
3 C diquark correlations should dominate the bound state contribution to the 
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three quark scattering matrix, and hence which should be retained in deriving 
and solving a Poincare covariant homogeneous Fadde'ev equation for baryons. 

The simple Goldstone-theorem-preserving rainbow-ladder kernel of Ref. |l6| 
can be used to illustrate this point. The model yields the following calculated 
diquark masses (isospin symmetry is assumed): 





(ud)o+ (us) Q + (uu)i+ (us)i+ (ss)i+ (uu)i- (us)i- (ss)i- 


m qq (GeV) 


0.74 0.88 0.95 1.05 1.13 1.47 1.53 1.64 



(2) 



and the results are relevant because the mass ordering is characteristic and 
model- independent, and lattice estimates, where available ]l7j, agree with the 
masses tabulated here. Equation (||) suggests that an accurate study of the 
nuclcon should retain the scalar and pseudovector correlations: (ud) +, (uu)i+, 
(ud)i+, (dd)i+, because for these diquarks m qq <m^, where is the nucleon 
mass, but may neglect other correlations. Furthermore, it is obvious from the 
angular momentum Clebsch-Gordon series: \ ® = ■§ and | ® 1 = | © §, that 
decuplet baryons are inaccessible without pseudovector diquark correlations. It 
is interesting to note that m( u[ ;) 0+ = 0.78 cf. 0.76 = tun I'm a and hence 
one might anticipate that the presence of diquark correlations in baryons can 
provide a straightforward explanation of the N-A mass-splitting and other like 
effects. These ideas were first enunciated in Refs. ]l8|,[l9) and Ref. |2(J provides 
a convincing demonstration of their efficacy. 

Explicit calculations; e.g., Ref. show that retaining only a scalar diquark 
correlation in the kernel of the nucleon's Fadde'ev equation provides insufficient 
binding to obtain the experimental nucleon mass: the best calculated value is 
typically ~ 40% too large. However, with the addition of a pseudovector diquark 
it is easy to simultaneously obtain [ ^2|pl| ] the experimental masses of the nucleon 
and A. Such calculations plainly verify the intuition that follows from simple 
mass-counting: the pseudovector diquarks are an important but subdominant 
clement of the nucleon's Fadde'ev amplitude (cf. the scalar diquark) whilst being 
the sole component of the A. 

The presence of diquark correlations in baryons also affects the predictions for 
scattering observables, which may therefore provide a means for experimentally 
verifying the ideas described above. For example, their presence provides a simple 
explanation of the neutron's nonzero electric form factor Jl2| : charge separation 
arising from a heavy (ud) diquark with electric charge | holding on to a relatively 
light, electric charge (— |) d-quark. And also a prediction for the ratio of the 
proton's valence-quark distributions: d/u := d v (x — > l)/u v (x — > 1), which can 
be measured in deep inelastic scattering |2^|. In this case, diquark correlations 
with differing masses in the nucleon's Fadde'ev amplitude are an immediate 
indication of the breaking of SU(6) symmetry, hence d/u ^ 1/2. Furthermore, 
if it were true that rri( qq ) jP 3 > TO( t id) 0+ > f° r au J P + > then d/u — 0. However, 
as we have seen, in reality the 1 + diquark is an important subdominant piece of 
the nucleon's Fadde'ev amplitude so that a realistic picture of diquarks in the 
nuclcon implies < d/u < i, with the actual value being a sensitive measure of 
the proton's pseudovector diquark fraction. 
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2 Superfluidity in Quark Matter 

We have outlined above the role and nature of diquark correlations in hadronic 
physics at zero temperature and density, and emphasised that diquarks are an 
idea as old as that of quarks themselves. Another phenomenon suggested imme- 
diately by the meson-diquark auxiliary-field effective action is that of diquark 
condensation; i.e., quark-quark Cooper pairing, which was first explored in this 
context using a simple version of the Nambu-Jona-Lasinio model . A chemi- 
cal potential promotes Cooper pairing in fermion systems and the possibility that 
such diquark pairing is exhibited in quark matter is also an old idea, early explo- 
rations of which employed (24] the rainbow-ladder truncation of the quark DSE 
(QCD gap equation). That interest in this possibility has been renewed is evident 
in a number of contributions to this volume p5| , p6|j2^ |. A quark-quark Cooper 
pair is a composite boson with both electric and colour charge, and hence super- 
fluidity in quark matter entails superconductivity and colour superconductivity. 
However, the last feature makes it difficult to identify an order parameter that 
can characterise a transition to the superfluid phase: the Cooper pair is gauge 
dependent and an order parameter is ideally describable by a gauge-invariant 
operator. This particularly inhibits an analysis of the phenomenon using lattice- 
QCD. 



2.1 Gap Equation 

Studies of the gap equation that suppress the possibility of diquark condensa- 
tion show that cold, sparse two-flavour QCD exhibits a nonzero quark-antiquark 
condensate: (qq) ^ 0. If it were otherwise then the 7r-meson would be almost 
as massive as the p-meson, which would yield a very different observable world. 
The quark condensate is undermined by increasing fx and T, and there is a large 
domain in the physical (upper-right) quadrant of the (/z, T)-plane for which 
(qq) = 0: for the purpose of exemplification, that domain can crudely be char- 
acterised as the set (see, e.g., Refs. p8pS|]3"oT ] ) : 

{(jj.,T) : (i 2 /ii 2 c +T 2 /T? > I, n,T > ; p, c - 0.3-0.4 GeV, T c ~ 0.15 GeV}. 

(3) 

Increasing temperature also opposes Cooper pairing. However, since increas- 
ing fi promotes it, there may be a (large-p,low-T)-subdomain in which quark 
matter exists in a superfluid phase. That domain, if it exists, is unlikely to be 
accessible at the Relativistic Heavy Ion Collider, because it operates in the high 
temperature regime, but may be realised in the core of compact astrophysical ob- 
jects, which could undergo a transition to superfluid quark matter as they cool. 



Possible signals accompanying such a transition are considered in Refs. [125 27:. 

It was observed in Ref. J3l| that a direct means of determining whether a 
SU C (N) gauge theory supports scalar diquark condensation is to study the gap 
equation satisfied by 
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Here T = 0, for illustrative simplicity and because temperature can only act to 
destabilise a condensate, and, with Uy = P4 + ifi, 

D(p, (i) = % 1 ■ p A(p 2 , Jfa ) + B(p 2 , u,^ ) + ^74 u [fl] C(p 2 , afa) ; (5) 

i.e., the inverse of the dressed-quark propagator in the absence of diquark pairing. 
(NB. For [i = 0, A, B and C are real functions.) It is one of the fundamental 
features of DSE studies that the existence of a nonzero quark condensate: (qq) ^ 
0, is signalled in the solution of the gap equation by B(p 2 ,u> 2 ^) ^00. 

In Eq. (|), {\\, i = l...n£, n£ = N C (N C - l)/2} are the antisymmetric 
generators of SU C (N C ) and C — 7274 is the charge conjugation matrix, 

C7*C+ = - 7;i ; [C, 75 ]=0, (6) 

where X t denotes the matrix transpose of X. The key new feature here is that 
diquark condensation is characterised by A l (p,/j,) ^ 0, for at least one i. That 
is clear if one considers the quark piece of the QCD Lagrangian density: L[q, q]. 
It is a scalar and hence L[q, q]* = L[q, q). Therefore L[q, q] oc L[q, q) + L[q, g]*, 
and it is a simple exercise to show that this sum, and hence the action, can be 
re-expressed in terms of a 2 x 2 diagonal matrix using the bispinor fields 

Q{x) := ( qix) .^2^) , Q{x) := (q(x) q(x) := q'Cr 2 ), (7) 

where {tj : i = 1,2, 3} are Pauli matrices that act on the isospin index It is 
plain upon inspection that a nonzero entry: d(x)"/s, in row-2-column-l of this 
action-matrix would act as a source for q t T 2 Cj$q; i.e., as a scalar diquark source. 

It is plain now that the explicit 2x2 matrix structure of T>(p, fi) in Eq. (|J) 
exhibits a quark bispinor index that is made with reference to Q{x), Q(x). This 
approach; i.e., employing a "matrix propagator" with "anomalous" off-diagonal 
elements, simply exploits the Gorkov-Nambu treatment of superconductivity in 
fermionic systems, which is explained in textbooks, e.g., Ref. |32|. It makes pos- 
sible a well-ordered treatment and makes unnecessary a truncated bosonisation, 
which in all but the simplest models is a procedure difficult to improve system- 
atically. 

The bispinor gap equation can be written in the form 

z>«,,,^„, P ,,, + ( Ti £ i ^ hiC£ ^ : ;;J, ct ), ( 8 , 

where the second term on the right-hand-side is just the bispinor self energy. 
Here, in the absence of a scalar diquark source term, 

T-^oiPi m) = (*7 ' V + m ) T Q — M t q > (9) 

with m the current-quark mass, and the additional Pauli matrices: {tq,o? = 
0, 1, 2, 3}, act on the bispinor indices. As we will see, the structure of £ij{p, y) 
specifies the theory and, in practice, also the approximation or truncation of it. 

1 We only consider theories with two light-flavours. Additional possibilities open if this 
restriction is lifted [ p5||2o| . 
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Two colour QCD (QC2D) provides an important and instructive example. In 
this case A l \\ = At 2 in Eq. (jl]), with \r c the generators of SU C (2), and it is 
useful to employ a modified bispinor 

Qi{x) := : S*q(x)) ' ( q{x) U x ) -IW^)- (10) 

Embedding the additional factor of r 2 in this way makes it possible to write the 
Lagrangian's fermion-gauge-boson interaction term as 



1 



Q2(x)- gl ^r^Q 2 (x)Al(x) (11) 

because SU C (2) is pseudoreal; i.e., t 2 (— r c ) = r c , and the fundamental and 
conjugate representations are equivalent; i.e., fermions and antifermions are prac- 
tically indistinguishable. (That the interaction term takes this form is easily seen 
using L[q,qf = L[q,q].) 

Using the pseudoreality of SU C (2) it can be shown that, for fi = and in the 
chiral limit, m = 0, the general solution of the bispinor gap equation is fjll 

V{p) = i"fpA{p 2 ) + V{-Tv)M(p 2 ), V(tt) = exp ji 75 ^ T l 7/j = V(-ir)- 1 , 

(12) 

where it 1 <---> 5 are arbitrary constants and 

{T 1,2,3 =t* q ® r f , T 4 = 7-1 ® r°, T 5 — Tq (E> r°} , {T 4 , T^"} = 2<5« , (13) 



S (P)= / A2( p 2)+ ^2( p2) : = -H-P^(p 2 ) + V(tt) ( t s (p 2 ). (14) 

[To illustrate this, note that inserting n = (0, 0, 0, 0, — j7r) produces an inverse 
bispinor propagator with the simple form in Eq. (0).] 

That the gap equation is satisfied for any constants ir e signals a vacuum 
degeneracy - it corresponds to a multidimensional "Mexican hat" structure of 
the theory's effective potential, as noted in a related context in Ref. (3j|. Con- 
sequently, if the interaction supports a mass gap, then that gap describes a 
five-parameter continuum of degenerate condensates: 

(Q 2 v( 7r )g 2 ) /o, (15) 

and there are 5 associated Goldstone bosons: 3 pions, a diquark and an anti- 
diquark. (Diquarks are the "baryons" of QC2D.) In the construction of Eq. ( p"2|) 
one has a simple elucidation of a necessary consequence of the Pauli-Giirsey 
symmetry of QC2D; i.e., the practical equivalence of particles and antiparticles. 
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Fig. 1. Illustration of the vertex-corrected gap equation, which is the next-to-leading- 
order in the systematic, symmetry-preserving truncation scheme of Ref. [Q. Retaining 
only the first two diagrams on the right-hand-side yields the dressed-rainbow trunca- 
tion. Each bispinor quark-gluon vertex is bare but the shaded circles mark quark and 
gluon 2-point functions that are dressed. The corresponding truncation in the relevant 
Bethe-Salpeter equations ensures the absence of diquark bound states in the strong 
interaction spectrum. (Adapted from Ref. pif .) 



For m ^ 0, the gap equation requires |u) trpQ [T l V] = 0, so that in this case 
only (Q2Q2) 7^ and now the spectrum contains five degenerate but massive 
pseudo-Goldstone bosons. This illustrates that a nonzero current-quark mass 
promotes a quark condensate and opposes diquark condensation. 

For /z ^ the general solution of the gap equation has the form 

In the absence of a diquark condensate; i.e., for A = 0, 

[U B (a),V(P,ti] - 0, U B (a) := e ia ^® T °f ; (17) 

i.e., baryon number is conserved in QC2D. This makes plain that the existence 
of a diquark condensate dynamically breaks this symmetry. 

To proceed we choose to be explicit and employ the dressed-rainbow trunca- 
tion of the gap equation, see Fig. |], with a model for the Landau gauge dressed- 
gluon propagator: 

kiiky \ r f 2 \ t- /; 2\ _ 64 4 „ 2 c-4/ 



9 Dfj,u{k) = U M „ - -^fj T 2 {k 2 ) , T*ip) = ff 8\k) . (18) 

This form for g 2 D fJ/l/ (k) was introduced |34| for the modelling of confinement 
in QCD but it is also appropriate here because the string tension in QC2D is 
nonzero, and that is represented implicitly in Eq. ( [l8| ) via the mass-scale r). 

Using Eq. @ we obtain an algebraic dressed-rainbow gap equation that, 
for p 2 — \p\ 2 + p 2 — 0, reads: 

A-l = ^rj 2 K {A{B* 2 -C* 2 y?) + A* \A\ 2 } , (19) 

(C-1),= I?K {C (B* 2 - C* 2 p 2 ) - C* \A\ 2 } , (20) 

B-m = if K {B {B* 2 - C* 2 li 2 ) + B* \A\ 2 } , (21) 
A =f) 2 K{A{\B\ 2 + \C\ 2 Li 2 ) + A\A\ 2 }, (22) 
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with K- 1 = \B 2 -C 2 fi 2 \ 2 + 2\A\ 2 (\B\ 2 + \C\ 2 n 2 ) + \A\ i . These equations possess 
a B <-> A symmetry when (m, /j,) — 0, which is a straightforward illustration of 
the vacuum degeneracy described above using the matrix V(tz). (Recall that for 
fi = 0, A, B and C are real functions.) They also exemplify the general result 
that A is real for all fi. Another exemplary result follows from a linearisation in 
H 2 : /i ^ acts to promote a nonzero value of A but oppose a nonzero value of 
B; i.e., a nonzero chemical potential plainly acts to promote Cooper pairing at 
the expense of (qq). 

For (m, fi) = the solution of the dressed-rainbow gap equation obtained 
using Eq. ( p~8|) is: 

^ 2) = c « ={j(i + yiTf), — , < 23 > 

X V) := B 2 (p 2 ) + A 2 (p 2 ) = h 2 - V, P 2 < f (24) 

[ 0, otherwise . 

As we have already mentioned, the dynamically generated mass function, _M(p 2 ), 
is tied to the existence of quark and/or diquark condensates, which can be il- 
lustrated by noting that (B = 0, A ^ 0) corresponds to ty = (0,0,0,0,^71") 
in Eq. fll5|); i.e., (Q2H5TQQ2) ^ 0, while (B 0,A = 0) corresponds to 
tt = (0,0,0,0,0); i.e., (Q 2 Q 2 ) ^ 0. 

The usual chiral, SUa{2), transformations are realised via 

V(p, fi) -> V(tt) X»(p, m) V(7r) , 7(tt) :^e t757r - T , tt = (tt 1 , tt 2 , tt 3 ) , (25) 

and therefore, since the anticommutator {T, T 4,5 } = 0, a diquark condensate 
does not dynamically break chiral symmetry. On the other hand, since [1, T] = 0, 
a quark condensate does dynamically break chiral symmetry. 

In addition, and of particular importance, is the feature that in combination 
with the momentum-dependent vector self energy the momentum-dependence of 
M.{p 2 ) ensures that the dressed-quark propagator does not have a Lehmann rep- 
resentation and hence can be interpreted as describing a confined quark pp|JlO| . 
The interplay between the scalar and vector self energies is the key to this realisa- 
tion of confinement. The qualitative features of this simple model's dressed-quark 
propagator have been confirmed in recent lattice-QCD simulations pq ] and the 
agreement between those simulations and more sophisticated DSE studies is 
semi-quantitative JljJ . 

In the steepest descent (or stationary phase) approximation the contribution 
of dressed-quarks to the thermodynamic pressure is 

Ps{^T) = jTrLn [PS- 1 ] - ^Tr [175]} , (26) 

where (3 — 1/T, and "Tr" and "Ln" are extensions of "tr" and "In" to matrix- 
valued functions. 
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The MIT Bag Model pictures the quarks in a baryon as occupying a spatial 
volume from which the nontrivial quark-condensed vacuum (scalar-field) has 
been expelled. Therefore, as observed in Refs. p3| ], the bag constant can be 
identified as the pressure difference between the (qq) / vacuum, the so-called 
Nambu-Goldstone phase in which chiral symmetry is dynamically broken, and 
the chirally symmetric no-condensate alternative, which is called the Wigner- 
Weyl vacuum. That difference is given by 

B B (p) ■= Pz{», S[B, A = 0])- psQi, S[B = 0, A = 0]) , (27) 

and it is, of course, /^-dependent because the vacuum evolves with changing /x. 
Bb also evolves with temperature and this necessary (p, Independence of the 
bag constant can have an important effect on quark star properties; e.g., reducing 
the maximum supportable mass of a quark matter star, as discussed in Ref. |3(|. 
If we define, by analogy, 

Ba(m) ■= PsQi, S[B = 0, A}) - Ps (fx, S[B = 0,A = 0}), (28) 

then the relative stability of the quark- and diquark-condensed phases is mea- 
sured by the pressure difference 

Sp(n) := B A (V) - Bb(v) ■ (29) 

For Sp(fi) > the diquark condensed phase is favoured. 
At (m = 0, p, = 0), 5p = 0, with 

Bb(O) =B A (0) = (0.092 f,) 4 . (30) 

This equality is a manifestation of the vacuum degeneracy identified above in 
connection with the matrix V(7r). However, 

with m = , dp > for all fi > , (31) 

which means that the Wigner-Weyl vacuum is unstable with respect to diquark 
condensation for all /i > and that the superfluid phase is favoured over 
the Nambu-Goldstone phase. 

Now, although the action for the /i ^ theory is invariant under 

Q 2 ^U B (a)Q 2 , Q 2 ^Q 2 U B (-a), (32) 

which is associated with baryon number conservation, the diquark condensate 
breaks this symmetry: 

(Q 2 H5TqQ 2 ) -> cos(2a) (Q 2 ij 5 TqQ 2 ) - sin(2a) {Q 2 ij 5 T~QQ 2 ) ; (33) 

i.e., it is a ground state that is not invariant under the transformation. Hence, 
for (m = 0, fi ^ 0), only one Goldstone mode remains. These symmetry breaking 
patterns and the concomitant numbers of Goldstone modes in QC2D are also 
described in Ref. fl3j| . 
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Fig. 2. Evolution of the critical chemical potential for diquark condensation as the 
current-quark mass is increased. The coordinate measures the magnitude of the current- 
quark mass through the mass of the theory's lightest excitation [a pseudo-Goldstone 
mode, as described after Eq. ([Hs])]. /i c and are measured in units of the model's 
mass scale, rj in Eq. (|li|): for m = 0, the vector meson mass is -j^V- 



For m ^ and small values of fj, the gap equation only admits a solution 
with A = 0; i.e., diquark condensation is blocked because the current-quark 
mass is a source of the quark condensate [see, e.g., Eq. ( f2l| ) and the comments 
after Eq. (^|)]. However, with increasing /i, the theory undergoes a transition to 
a phase in which the diquark condensate is nonzero. We identify the transition 
as second order because the diquark condensate falls continuously to zero as 
fj, — > /i+ , where fi c is the critical chemical potential. In Fig. ^| we plot the critical 
chemical potential as a function of /2, where, to sidestep solving the Bethe- 
Salpeter equation, was obtained using a Gell-Mann-Oakes-Renner-like mass 
formula, Eqs. (16)-(18) in Ref. p9[ , which follows [B8| from the axial-vector 
Ward-Takahashi identity. From the figure it is clear that this simple model of 
QC2D exhibits the relation 

Mc = 2 m ir i (34) 

which is anticipated for QCD-like theories with pseudoreal fermions |${|. We note 
that the deviation from Eq. ([54]) at larger values of results from neglecting 
0(r?7 2 )-corrections in the mass formula. This omission leads to an underestimate 
of the pion mass |4C|] , which is responsible for the upward deflection of the 
calculated results evident in Fig. g. 

In exemplifying these features we have employed the rainbow-ladder trunca- 
tion. However, improving on that will only yield quantitative changes of < 20% 
in the results because the pseudoreality of QC2D and the equal dimension of the 
colour and bispinor spaces, which underly the theory's Pauli-Gurscy symmetry, 
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ensure that the entire discussion remains qualitatively unchanged. In particular, 
the results of Fig. || and Eq. (|34|), being tied to chiral symmetry, remain un- 
changed because at least one systematic, chiral symmetry preserving truncation 
scheme exists Q . 



4 Three Colours 



The exploration of superfluidity in true QCD encounters two differences: the 
dimension of the colour space is greater than that of the bispinor space and the 
fundamental and conjugate representations of the gauge group are not equiva- 
lent. The latter is of obvious importance because it entails that the quark-quark 
and quark- antiquark scattering matrices are qualitatively different. 

n£ — 3 in QCD and hence in canvassing superfluidity it is necessary to choose 
a direction for the condensate in colour space;f] e.g., A l \\ = A A 2 in Eq. (|]), so 
that 

V(p, fx) 



D n (p,f,)P n +D ± (p,fi)P ± 


A( P ,nh 5 \ 2 


-A(p, -/x) 75 A 2 


CD\\ (-p, M)CtP|, + CD ± (p, n)&P± . 



where P|| = (A 2 ) 2 , P± + P\\ — diag(l, 1,1), and D\\, D± are defined via obvious 
generalisations of Eqs. (||), (||). In Eq. (|35| ) the evident, demarcated block struc- 
ture makes explicit the bispinor index: each block is a 3 x 3 colour matrix and 
the subscripts: ||, _L, indicate whether or not the subspace is accessible via A2. 

The bispinors associated with this representation are given in Eqs. (Q) and 
in this case the Lagrangian's quark-gluon interaction term is 



Q(x)igr«Q(x)A«(x), 1% 














(36) 



It is instructive to compare this with Eq. ([ll] ) : with three colours the interaction 
term is not proportional to the identity matrix in the bispinor space, Tq. This 
makes plain the inequivalence of the fundamental and conjugate fermion repre- 
sentations of SU C (3), which entails that quark- antiquark scattering is different 
from quark-quark scattering. 

It is straightforward to derive the gap equation at arbitrary order in the 
truncation scheme of Ref. [[lJJ and it is important to note that because 

D\\(p,n)P\\+D x (p,iJ l )P ± = (37) 
A { §D|| (p, n) + §D ± (p, m)} + 75A 8 {D|| (p, M ) - D ± (p, p) } 

the interaction: r^S(p, fi)T^ , necessarily couples the ||- and _L-components. Ref- 
erence plf explored the possibility of diquark condensation in QCD using both 



2 It is this selection of a direction in colour space that opens the possibility for 
colour-flavour locked diquark condensation in a theory with three effectively-massless 



quarks; i.e., current-quark masses p^ , pr| . 
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the rainbow and vertex-corrected gap equation, illustrated in Fig. El with 

g 2 D^(k) = (s^- T{k 2 ) , Hk 2 ) = 4^ 4 ry 2 5\k) . (38) 



For (m, ji) — the rainbow-ladder truncation yields 



rn 



2 _ 1 „2 
P 2 



i^ 2 , (qq)° = (0.11 t?) 3 , B B (fi = 0) = (0.10 t?) 4 , (39) 



and momentum-dependent vector self energies, which lead to an interaction be- 
tween the ||- and _L-components of T> that blocks diquark condensation. This is 
in spite of the fact that A a A 2 (— A a )* = §A a A a , which entails that the rainbow- 
truncation quark-quark scattering kernel is purely attractive and strong enough 
to produce diquark bound states jl3| . (Remember that in the colour singlet me- 
son channel the rainbow-ladder truncation gives the colour coupling A a A a ; i.e., 
an interaction with the same sign but twice as strong.) 

For fi and in the absence of diquark condensation this model and trun- 
cation exhibits |^8| coincident, first order chiral symmetry restoring and decon- 
fining transitions at 

M c B rai;bow = 0.28 ry = 0.3 GeV, (40) 



with r\ = 1.06 GeV fixed by fitting the ra/0 vector meson mass |14 . 

For (to = 0, /i =/= 0), however, the rainbow-truncation gap equation admits 
a solution with A(p, fi) ^ and B{p,n) = 0. The pressure difference, 5p([i) in 
Eq. (|29|), is again the way to determine whether the stable ground state is the 
Nambu-Goldstone or superfluid phase. With increasing fi, Bb(^) decreases, very 
slowly at first, and Ba(p-) increases rapidly from zero. That evolution continues 
until 

Vc, Ta°nbow = °' 25 V = °' 89 ^'rftabow . ( 41 ) 

where Ba(^) becomes greater-than Bb{^)- This signals a first order transition 
to the superfluid ground state and at the boundary 

(Qi l5 T 2 X 2 Q)^B : o,^ = (0.65) 3 (QQ)^ . (42) 

Since Ba {fJ-) > for all /z > there is no intermediate domain of \x in which all 
condensates vanish. 

The solution of the rainbow gap equation: A(p,nf'^), which is real and 
characterises the diquark gap, is plotted in Fig. |j| It vanishes at p 2 — as a 

consequence of the 1 1 — I coupling that blocked diquark condensation at fx = 0, 

and also at large p 2 , which is a manifestation of this simple model's version of 
asymptotic freedom. 

The chemical potential^] at which the switch to the superfluid ground state 
occurs, Eq. (f4l|), is consistent with other estimates made using models compa- 

3 We note that in a two flavour free-quark gas at fx = 0.3 GeV the baryon number 
density is 1.5 po, where po = 0.16 fm -3 . In the same system at fi = 0.55 GeV the 
baryon number density is > 10 po. For comparison, the central core density expected 
in a 1.4Mq neutron star is 3.6-4.1 po [^lj. Arguments valid at "asymptotically large" 
quark chemical potential are therefore unlikely to be relevant to experimentally or 
observationally accessible systems. 
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Fig. 3. Dashed line: A(z, ^c' A ) obtained in rainbow truncation with the QCD model 
defined via Eq. d3g|), plotted for a = as a function of p, where z = p (0, 0, sin a, i/i + 
cos a). As fi increases, the peak position shifts to larger values of p and the peak height 
increases. Solid line: A(z,fi = 0) obtained as the solution of Eq. (fl3|), the vertex- 
corrected gap equation, also with a = 0. (Adapted from Ref. 



rable to the rainbow-truncation class 1 25 2r^ , ^2] , ^3|j44| , as is the large magnitude 
of the gap at this point |^||,||,||. 

A question that now arises is: How sensitive is this phenomenon to the nature 
of the quark-quark interaction? As we discussed in connection with Eq. (^J), the 
inhomogeneous dressed-ladder BSE exhibits particle-like singularities in the + 
diquark channels and such states do not exist in the strong interaction spectrum. 
Does diquark condensation persist when a truncation of the gap equation is 
employed that does not correspond to a BSE whose solutions exhibit diquark 
bound states? The vertex corrected gap equation, 



P(p,M) = (p. M) 



(43) 



J^ 4 r; s(p, m) rl s(p, r; s( P , M ) r b a . 



which is depicted in Fig. [3, is just such a truncation and it was also studied in 
Ref. g. 

In this case there is a A ^ solution even for /i = 0, which is illustrated in 
Fig. ||, and using the interaction of Eq. ( p^ ) 

ml = (1.1) m 2 p laddcr , (QQ) = (1.0) 3 (QQ) rainbow , B B = (l.l) 4 B^ ainbow , (44) 

where the rainbow-ladder results are given in Eqs. (|39|), and 

(QiTo T Q^ 2 Q) = (0-48) 3 (QQ) , B A = (0.42) 4 B B . (45) 
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The last result shows, unsurprisingly, that the Nambu-Goldstone phase is fa- 
voured at fj, = 0. Precluding diquark condensation, the solution of the vertex- 
corrected gap equation exhibits coincident, first order chiral symmetry restoring 
and deconfincmcnt transitions at 

B,A=0 n nn B,A=Q fA „s 

H c ' =0.77/^ c ; rainbow . (46) 

Admitting diquark condensation, however, the /Li-dependence of the bag con- 
stants again shows there is a first order transition to the superfiuid phase, here 
at 

tf=o,A = .63 M f 4 =°, with (Qi l5 Tl\ 2 Q)^ 63 ^,A= = (0.51) 3 (QQ) M=0 . 

(47) 

(NB. This discussion is still for m = 0. We saw at the end of Sec. || what effects 
to anticipate at m / 0.) Thus the material step of employing a truncation that 
eliminates diquark bound states leads only to small quantitative changes in the 
quantities characterising the still extant superfiuid phase; e.g., reductions in the 
magnitude of both the critical chemical potential for the transition to superfiuid 
quark matter and the gap. Hence scalar diquark condensation appears to be 
a robust phenomenon. One caveat to bear in mind, however, is that the gap 
equation studies conducted hitherto do not obviate the question of whether the 
diquark condensed phase is stable with-respect-to dinucleon condensation f45|| , 
which requires further attention. 

Heating causes the diquark condensate to evaporate. Existing studies suggest 

that it will disappear for T > 60-100 MeV HI|,||,|4). However ' such tem P er " 
atures are high relative to that anticipated inside dense astrophysical objects, 
which may indeed therefore provide an environment for detecting quark matter 
superfluidity. 



5 Summary 

The idea that diquark correlations play an important role in strong interaction 
physics is an old one. However, modern computational resources and theoretical 
techniques make possible a more thorough and quantitative exploration of the 
merits of this idea and its realisation in QCD. These advances are in part respon- 
sible for the contemporary resurgence of interest in all aspects of diquark-related 
phenomena. 

Herein we have attempted to provide a qualitative understanding of the na- 
ture of diquark condensation using exemplary, algebraic models, and focusing 
on two flavour theories for simplicity. 

The gap equation is a primary tool in all studies of pairing. Using the special 
case of 2-colour QCD, QC2D, we illustrated via an analysis of the gap equation 
how a nonzero chemical potential promotes Cooper pairing and how that pairing 
can overwhelm a source for quark-antiquark condensation, such as a fermion 
current-mass. As we saw, the pseudoreality of SU(N C = 2) entails that QC2D 
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has a number of special symmetry properties, which dramatically affect the 
spectrum. 

Turning to QCD itself, we saw that one can expect a nonzero quark conden- 
sate at zero chemical potential: a := — {qq) 1 ^ 3 / 0, to give way to a diquark 
condensate when the chemical potential exceeds ~ 2a, and at this point the 
diquark gap is ~ a/2. The diquark condensate melts when the temperature ex- 
ceeds ~ 60 - 100 MeV; i.e., one-third to one-half of the chiral symmetry restoring 
temperature in two-flavour QCD. These features are model-independent in the 
sense that the many, disparate models applied recently to the problem yield 
results in semi-quantitative agreement. 
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